High-risk and silent clonal hematopoietic genotypes in patients with nonhematologic cancer.

ABSTRACT: Clonal hematopoiesis (CH) identified by somatic gene variants with variant allele fraction (VAF) ≥ 2% is associated with an increased risk of hematologic malignancy. However, CH defined by a broader set of genotypes and lower VAFs is ubiquitous in older individuals. To improve our understanding of the relationship between CH genotype and risk of hematologic malignancy, we analyzed data from 42 714 patients who underwent blood sequencing as a normal comparator for nonhematologic tumor testing using a large cancer-related gene panel. We cataloged hematologic malignancies in this cohort using natural language processing and manual curation of medical records. We found that some CH genotypes including JAK2, RUNX1, and XPO1 variants were associated with high hematologic malignancy risk. Chronic disease was predicted better than acute disease suggesting the influence of length bias. To better understand the implications of hematopoietic clonality independent of mutational function, we evaluated a set of silent synonymous and noncoding mutations. We found that silent CH, particularly when multiple variants were present or VAF was high, was associated with increased risk of hematologic malignancy. We tracked expansion of CH mutations in 26 hematologic malignancies sequenced with the same platform. JAK2 and TP53 VAF consistently expanded at disease onset, whereas DNMT3A and silent CH VAFs mostly decreased. These data inform the clinical and biological interpretation of CH in the context of nonhematologic cancer.

Associations Between Cancer Predisposition Mutations and Clonal Hematopoiesis in Patients With Solid Tumors.

PURPOSE: Clonal hematopoiesis (CH), the expansion of clones in the hematopoietic system, has been linked to different internal and external features such as aging, genetic ancestry, smoking, and oncologic treatment. However, the interplay between mutations in known cancer predisposition genes and CH has not been thoroughly examined in patients with solid tumors. METHODS: We used prospective tumor-blood paired sequencing data from 46,906 patients who underwent Memorial Sloan Kettering-Integrated Mutation Profiling of Actionable Cancer Targets (MSK-IMPACT) testing to interrogate the associations between CH and rare pathogenic or likely pathogenic (P/LP) germline variants. RESULTS: We observed an enrichment of CH-positive patients among those carrying P/LP germline mutations and identified a significant association between P/LP germline variants in ATM and CH. Germline and CH comutation patterns in ATM, TP53, and CHEK2 suggested biallelic inactivation as a potential mediator of clonal expansion. Moreover, we observed that CH-PPM1D mutations, similar to somatic tumor-associated PPM1D mutations, were depleted in patients with P/LP germline mutations in the DNA damage response (DDR) genes ATM, CHEK2, and TP53. Patients with solid tumors and harboring P/LP germline mutations, CH mutations, and mosaicism chromosomal alterations might be at an increased risk of developing secondary leukemia while germline variants in TP53 were identified as an independent risk factor (hazard ratio, 36; P < .001) for secondary leukemias. CONCLUSION: Our results suggest a close relationship between inherited variants and CH mutations within the DDR genes in patients with solid tumors. Associations identified in this study might translate into enhanced clinical surveillance for CH and associated comorbidities in patients with cancer harboring these germline mutations.

Germline mutations in WNK2 could be associated with serrated polyposis syndrome.

BACKGROUND: Patients with serrated polyposis syndrome (SPS) have multiple and/or large serrated colonic polyps and higher risk for colorectal cancer. SPS inherited genetic basis is mostly unknown. We aimed to identify new germline predisposition factors for SPS by functionally evaluating a candidate gene and replicating it in additional SPS cohorts. METHODS: After a previous whole-exome sequencing in 39 SPS patients from 16 families (discovery cohort), we sequenced specific genes in an independent validation cohort of 211 unrelated SPS cases. Additional external replication was also available in 297 SPS cases. The WNK2 gene was disrupted in HT-29 cells by gene editing, and WNK2 variants were transfected using a lentiviral delivery system. Cells were analysed by immunoblots, real-time PCR and functional assays monitoring the mitogen-activated protein kinase (MAPK) pathway, cell cycle progression, survival and adhesion. RESULTS: We identified 2 rare germline variants in the WNK2 gene in the discovery cohort, 3 additional variants in the validation cohort and 10 other variants in the external cohorts. Variants c.2105C>T (p.Pro702Leu), c.4820C>T (p.Ala1607Val) and c.6157G>A (p.Val2053Ile) were functionally characterised, displaying higher levels of phospho-PAK1/2, phospho-ERK1/2, CCND1, clonogenic capacity and MMP2. CONCLUSION: After whole-exome sequencing in SPS cases with familial aggregation and replication of results in additional cohorts, we identified rare germline variants in the WNK2 gene. Functional studies suggested germline WNK2 variants affect protein function in the context of the MAPK pathway, a molecular hallmark in this disease.

PCIG: a web-based application to explore immune-genomics interactions across cancer types.

MOTIVATION: Genomic alterations can modulate the tumor immunophenotype depending on their nature and tissue of origin. Although this immune-genomic interaction may shape disease progression and response to immunotherapy, the factors governing such dynamics and the influence of each tissue-specific context remain poorly understood. RESULTS: Here, we have developed the PanCancer ImmunoGenomics (PCIG) tool, a web-based resource that provides researchers with the opportunity to mine immunome-genome relationships across several cancer types using data from the Pan-Cancer Analysis of Whole-Genomes (PCAWG) study, which comprises >2,600 samples spanning across 20 different cancer primary sites. PCIG yields an integrative analysis of the crosstalk between somatic genomic alterations and different immune features, thus helping to understand immune response-related processes. AVAILABILITY AND IMPLEMENTATION: PCIG is freely available at https://pcig.vhio.net and is supported by all major web browsers. PCIG was developed with Django, which is a Python-based free and open-source framework, and it uses SQL Server as a relational database management system. The code is freely available for download at GitHub https://github.com/AnnaPG/PCIG and in its online supplementary material. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Integrated multi-omics profiling of high-grade estrogen receptor-positive, HER2-negative breast cancer.

Estrogen receptor-positive and human epidermal growth factor receptor 2-negative (ER+ HER2- ) breast cancer accounts for ~ 60-70% of all cases of invasive breast carcinoma. High-grade ER+ HER2- tumors respond poorly to endocrine therapy. In this study, we systematically analyzed clinical and multi-omics data to find potential strategies for personalized therapy of patients with high-grade ER+ HER2- disease. Six different cohorts were analyzed, for which multi-omics data were available. Grade III ER+ HER2- cases harbored higher proportions of large tumor size (> 5 cm), lymph node metastasis, chemotherapy use, and luminal B subtypes defined by PAM50, as compared with grade I/II tumors. DNA methylation (HM450) data and methylation-specific PCR indicated that the cg18629132 locus in the MKI67 promoter was hypermethylated in grade I/II cases and normal tissue, but hypomethylated in grade III cases or triple-negative breast cancer, resulting in higher expression of MKI67. Mutations in ESR1 and TP53 were detected in post-endocrine treatment metastatic samples at a higher rate than in treatment-naive tumors in grade III cases. We identified 42 and 20 focal copy number events in nonmetastatic and metastatic high-grade ER+ HER2- cases, respectively, with either MYC or MDM2 amplification representing an independent prognostic event in grade III cases. Transcriptional profiling within grade III tumors highlighted ER signaling downregulation and upregulation of immune-related pathways in non-luminal-like tumors defined by PAM50. Recursive partitioning analysis was employed to construct a decision tree of an endocrine-resistant subgroup (GATA3-negative and AGR-negative) of two genes that was validated by immunohistochemistry in a Chinese cohort. All together, these data suggest that grade III ER+ HER2- tumors have distinct clinical and molecular characteristics compared with low-grade tumors, particularly in cases with non-luminal-like biology. Due to the dismal prognosis in this group, clinical trials are warranted to test the efficacy of potential novel therapies.

Comprehensive Genomic Characterization of Fifteen Early-Onset Lynch-Like Syndrome Colorectal Cancers.

Lynch-like syndrome (LLS) is an increasingly common clinical challenge with an underlying molecular basis mostly unknown. To shed light onto it, we focused on a very young LLS early-onset colorectal cancer (CRC) cohort (diagnosis ≤ 40 y.o.), performing germline and tumor whole-exome sequencing (WES) of 15 patients, and additionally analyzing their corresponding tumor mutational burden (TMB) and mutational signatures. We identified four cases (27%) with double somatic putative variants in mismatch repair (MMR) core genes, as well as three additional cases (20%) with double MSH3 somatic alterations in tumors with unexplained MSH2/MSH6 loss of expression, and two cases (13%) with POLD1 potential biallelic alterations. Average TMB was significantly higher for LLS cases with double somatic alterations. Lastly, nine predicted deleterious variants in genes involved in the DNA repair functions and/or previously associated with CRC were found in nine probands, four of which also showed MMR biallelic somatic inactivation. In conclusion, we contribute new insights into LLS CRC, postulating MSH3 and POLD1 double somatic alterations as an underlying cause of a microsatellite instability (MSI) phenotype, proposing intrinsic biological differences between LLS with and without somatic alterations, and suggesting new predisposing candidate genes in this scenario.

Germline and Somatic Whole-Exome Sequencing Identifies New Candidate Genes Involved in Familial Predisposition to Serrated Polyposis Syndrome.

The serrated polyposis syndrome (SPS) is the most common and yet underdiagnosed colorectal polyposis syndrome. It is characterized by multiple and/or large colonic serrated polyps and a higher associated risk for colorectal cancer (CRC). The main objective of this study was to identify new candidate genes involved in the germline predisposition to SPS/CRC. Thirty-nine SPS patients from 16 families (≥2 patients per family) were recruited without alterations in well-known hereditary CRC genes, and germline and somatic whole-exome sequencing were performed. Germline rare variants with plausible pathogenicity, located in genes involved in cancer development, senescence and epigenetic regulation were selected. Somatic mutational profiling and signature analysis was pursued in one sample per family, when possible. After data filtering, ANXA10, ASXL1, CFTR, DOT1L, HIC1, INO80, KLF3, MCM3AP, MCM8, PDLIM2, POLD1, TP53BP1, WNK2 and WRN were highlighted as the more promising candidate genes for SPS germline predisposition with potentially pathogenic variants shared within families. Somatic analysis characterized mutational profiles in advanced serrated polyps/tumors, revealing a high proportion of hypermutated samples, with a prevalence of clock-like mutational signatures in most samples and the presence of DNA mismatch repair-defective signatures in some cases. In conclusion, we identified new candidate genes to be involved in familial SPS. Further functional studies and replication in additional cohorts are required to confirm the selected candidates.

Interplay between chromosomal alterations and gene mutations shapes the evolutionary trajectory of clonal hematopoiesis.

Stably acquired mutations in hematopoietic cells represent substrates of selection that may lead to clonal hematopoiesis (CH), a common state in cancer patients that is associated with a heightened risk of leukemia development. Owing to technical and sample size limitations, most CH studies have characterized gene mutations or mosaic chromosomal alterations (mCAs) individually. Here we leverage peripheral blood sequencing data from 32,442 cancer patients to jointly characterize gene mutations (n = 14,789) and mCAs (n = 383) in CH. Recurrent composite genotypes resembling known genetic interactions in leukemia genomes underlie 23% of all detected autosomal alterations, indicating that these selection mechanisms are operative early in clonal evolution. CH with composite genotypes defines a patient group at high risk of leukemia progression (3-year cumulative incidence 14.6%, CI: 7-22%). Multivariable analysis identifies mCA as an independent risk factor for leukemia development (HR = 14, 95% CI: 6-33, P < 0.001). Our results suggest that mCA should be considered in conjunction with gene mutations in the surveillance of patients at risk of hematologic neoplasms.

Copy-Number Alteration Burden Differentially Impacts Immune Profiles and Molecular Features of Hepatocellular Carcinoma.

PURPOSE: Chromosomal instability is a hallmark of cancer that results in broad and focal copy-number alterations (CNAs), two events associated with distinct molecular, immunologic, and clinical features. In hepatocellular carcinoma (HCC), the role of CNAs has not been thoroughly assessed. Thus, we dissected the impact of CNA burdens on HCC molecular and immune features. EXPERIMENTAL DESIGN: We analyzed SNP array data from 452 paired tumor/adjacent resected HCCs and 25 dysplastic nodules. For each sample, broad and focal CNA burdens were quantified using CNApp, and the resulting broad scores (BS) and focal scores (FS) were correlated with transcriptomic, mutational, and methylation profiles, tumor immune composition, and clinicopathologic data. RESULTS: HCCs with low broad CNA burdens (defined as BS ≤ 4; 17%) presented high inflammation, active infiltrate signaling, high cytolytic activity, and enrichment of the "HCC immune class" and gene signatures related to antigen presentation. Conversely, tumors with chromosomal instability (high broad CNA loads, BS ≥ 11; 40%), displayed immune-excluded traits and were linked to proliferation, TP53 dysfunction, and DNA repair. Candidate determinants of the low cytotoxicity and immune exclusion features of high-BS tumors included alterations in antigen-presenting machinery (i.e., HLA), widespread hypomethylation, and decreased rates of observed/expected neoantigenic mutations. High FSs were independent of tumor immune features, but were related to proliferation, TP53 dysfunction, and progenitor cell traits. CONCLUSIONS: HCCs with high chromosomal instability exhibit features of immune exclusion, whereas tumors displaying low burdens of broad CNAs present an immune active profile. These CNA scores can be tested to predict response to immunotherapies.

Germline Mutations in FAF1 Are Associated With Hereditary Colorectal Cancer.

BACKGROUND & AIMS: A significant proportion of colorectal cancer (CRC) cases have familial aggregation but little is known about the genetic factors that contribute to these cases. We performed an exhaustive functional characterization of genetic variants associated with familial CRC. METHODS: We performed whole-exome sequencing analyses of 75 patients from 40 families with a history of CRC (including early-onset cases) of an unknown germline basis (discovery cohort). We also sequenced specific genes in DNA from an external replication cohort of 473 families, including 488 patients with colorectal tumors that had normal expression of mismatch repair proteins (validation cohort). We disrupted the Fas-associated factor 1 gene (FAF1) in DLD-1 CRC cells using CRISPR/Cas9 gene editing; some cells were transfected with plasmids that express FAF1 missense variants. Cells were analyzed by immunoblots, quantitative real-time polymerase chain reaction, and functional assays monitoring apoptosis, proliferation, and assays for Wnt signaling or nuclear factor (NF)-kappa-B activity. RESULTS: We identified predicted pathogenic variant in the FAF1 gene (c.1111G>A; p.Asp371Asn) in the discovery cohort; it was present in 4 patients of the same family. We identified a second variant in FAF1 in the validation cohort (c.254G>C; p.Arg85Pro). Both variants encoded unstable FAF1 proteins. Expression of these variants in CRC cells caused them to become resistant to apoptosis, accumulate beta-catenin in the cytoplasm, and translocate NF-kappa-B to the nucleus. CONCLUSIONS: In whole-exome sequencing analyses of patients from families with a history of CRC, we identified variants in FAF1 that associate with development of CRC. These variants encode unstable forms of FAF1 that increase resistance of CRC cells to apoptosis and increase activity of beta-catenin and NF-kappa-B.

CNApp, a tool for the quantification of copy number alterations and integrative analysis revealing clinical implications.

Somatic copy number alterations (CNAs) are a hallmark of cancer, but their role in tumorigenesis and clinical relevance remain largely unclear. Here, we developed CNApp, a web-based tool that allows a comprehensive exploration of CNAs by using purity-corrected segmented data from multiple genomic platforms. CNApp generates genome-wide profiles, computes CNA scores for broad, focal and global CNA burdens, and uses machine learning-based predictions to classify samples. We applied CNApp to the TCGA pan-cancer dataset of 10,635 genomes showing that CNAs classify cancer types according to their tissue-of-origin, and that each cancer type shows specific ranges of broad and focal CNA scores. Moreover, CNApp reproduces recurrent CNAs in hepatocellular carcinoma and predicts colon cancer molecular subtypes and microsatellite instability based on broad CNA scores and discrete genomic imbalances. In summary, CNApp facilitates CNA-driven research by providing a unique framework to identify relevant clinical implications. CNApp is hosted at https://tools.idibaps.org/CNApp/.

In most cases, human cells contain two copies of each of their genes, yet sometimes this can change, an effect called copy number alteration (CNA). Cancer is a genetic disease and thus, studying the DNA from tumor samples is crucial to improving diagnosis and choosing the right treatment. Most tumors contain cells with CNAs; however, the impact of CNAs in cancer progression is poorly understood. CNAs can be studied by examining the genome of tumor cells and finding which regions display an unusual number of copies. It may also be possible to gather information about different cancer types by analyzing the CNAs in a tumor, but this approach requires the analysis of large amounts of data. To aid the analysis of CNAs in cancer cells, Franch-Expósito, Bassaganyas et al. have created an online tool called CNApp, which is able to identify and count CNAs in genomic data and link them to features associated with different cancers. The hope is that a better understanding of the effect of CNAs in cancer could help better diagnose cancers, and improve outcomes for patients. Potentially, this could also predict what type of treatment would work better for a specific tumor. Besides, by using a machine-learning approach, the tool can also make predictions about specific cancer subtypes in order to facilitate clinical decisions. Franch-Expósito, Bassaganyas et al. tested CNApp using previously existing cancer data from 33 different cancer types to show how CNApp can help the interpretation of CNAs in cancer. Moreover, CNApp can also use CNAs to identify different types of bowel (colorectal) cancer in a way that could help doctors to make decisions about treatment. Together these findings show that CNApp provides an adaptable and accessible research tool for the study of cancer genomics, which could provide opportunities to inform medical procedures.

Using linkage studies combined with whole-exome sequencing to identify novel candidate genes for familial colorectal cancer.

Colorectal cancer (CRC) is a complex disorder for which the majority of the underlying germline predisposition factors remain still unidentified. Here, we combined whole-exome sequencing (WES) and linkage analysis in families with multiple relatives affected by CRC to identify candidate genes harboring rare variants with potential high-penetrance effects. Forty-seven affected subjects from 18 extended CRC families underwent WES. Genome-wide linkage analysis was performed under linear and exponential models. Suggestive linkage peaks were identified on chromosomes 1q22-q24.2 (maxSNP = rs2134095; LODlinear = 2.38, LODexp = 2.196), 7q31.2-q34 (maxSNP = rs6953296; LODlinear = 2.197, LODexp = 2.149) and 10q21.2-q23.1 (maxSNP = rs1904589; LODlinear = 1.445, LODexp = 2.195). These linkage signals were replicated in 10 independent sets of random markers from each of these regions. To assess the contribution of rare variants predicted to be pathogenic, we performed a family-based segregation test with 89 rare variants predicted to be deleterious from 78 genes under the linkage intervals. This analysis showed significant segregation of rare variants with CRC in 18 genes (weighted p-value > 0.0028). Protein network analysis and functional evaluation were used to suggest a plausible candidate gene for germline CRC predisposition. Etiologic rare variants implicated in cancer germline predisposition may be identified by combining traditional linkage with WES data. This approach can be used with already available NGS data from families with several sequenced members to further identify candidate genes involved germline predisposition to disease. This approach resulted in one candidate gene associated with increased risk of CRC but needs evidence from further studies.

Identification of a Novel Candidate Gene for Serrated Polyposis Syndrome Germline Predisposition by Performing Linkage Analysis Combined With Whole-Exome Sequencing.

OBJECTIVES: Serrated polyposis syndrome (SPS) is a complex disorder with a high risk of colorectal cancer for which the germline factors remain largely unknown. Here, we combined whole-exome sequencing (WES) and linkage studies in families with multiple members affected by SPS to identify candidate genes harboring rare variants with higher penetrance effects. METHODS: Thirty-nine affected subjects from 16 extended SPS families underwent WES. Genome-wide linkage analysis was performed under linear and exponential models. The contribution of rare coding variants selected to be highly pathogenic was assessed using the gene-based segregation test. RESULTS: A significant linkage peak was identified on chromosome 3p25.2-p22.3 (maxSNP = rs2293787; LODlinear = 2.311, LODexp = 2.11), which logarithm of the odds (LOD) score increased after fine mapping for the same marker (maxSNP = rs2293787; LODlinear = 2.4, LODexp = 2.25). This linkage signal was replicated in 10 independent sets of random markers from this locus. To assess the contribution of rare variants predicted to be pathogenic, we performed a family-based segregation test with 11 rare variants predicted to be deleterious from 10 genes under the linkage intervals. This analysis showed significant segregation of rare variants with SPS in CAPT7, TMEM43, NGLY1, and FBLN2 genes (weighted P value > 0.007). DISCUSSION: Protein network analysis suggested FBLN2 as the most plausible candidate genes for germline SPS predisposition. Etiologic rare variants implicated in disease predisposition may be identified by combining traditional linkage with WES data. This powerful approach was effective for the identification of a new candidate gene for hereditary SPS.

The landscape of genomic copy number alterations in colorectal cancer and their consequences on gene expression levels and disease outcome.

Aneuploidy, the unbalanced state of the chromosome content, represents a hallmark of most solid tumors, including colorectal cancer. Such aneuploidies result in tumor specific genomic imbalances, which emerge in premalignant precursor lesions. Moreover, increasing levels of chromosomal instability have been observed in adenocarcinomas and are maintained in distant metastases. A number of studies have systematically integrated copy number alterations with gene expression changes in primary carcinomas, cell lines, and experimental models of aneuploidy. In fact, chromosomal aneuploidies target a number of genes conferring a selective advantage for the metabolism of the cancer cell. Copy number alterations not only have a positive correlation with expression changes of the majority of genes on the altered genomic segment, but also have effects on the transcriptional levels of genes genome-wide. Finally, copy number alterations have been associated with disease outcome; nevertheless, the translational applicability in clinical practice requires further studies. Here, we (i) review the spectrum of genetic alterations that lead to colorectal cancer, (ii) describe the most frequent copy number alterations at different stages of colorectal carcinogenesis, (iii) exemplify their positive correlation with gene expression levels, and (iv) discuss copy number alterations that are potentially involved in disease outcome of individual patients.

Somatic mutational signatures in polyposis and colorectal cancer.

The somatic mutation spectrum imprinted in the genome of a tumor represents the mutational processes that have been active in that tumor. Large sequencing efforts in various cancer types have resulted in the identification of multiple mutational signatures, of which several have been linked to specific biological mechanisms. Several pan-cancer mutational signatures have been identified, while other signatures are only found in specific tissue types. Research on tumors from individuals with specific DNA repair defects has led to links between specific mutational signatures and mutational processes. Studying mutational signatures in cancers that are likely the result of a genetic predisposition may represent an interesting strategy to identify constitutional DNA repair defects, including those underlying polyposis and colorectal cancer.

Approaches to functionally validate candidate genetic variants involved in colorectal cancer predisposition.

Most next generation sequencing (NGS) studies identified candidate genetic variants predisposing to colorectal cancer (CRC) but do not tackle its functional interpretation to unequivocally recognize a new hereditary CRC gene. Besides, germline variants in already established hereditary CRC-predisposing genes or somatic variants share the same need when trying to categorize those with relevant significance. Functional genomics approaches have an important role in identifying the causal links between genetic architecture and phenotypes, in order to decipher cellular function in health and disease. Therefore, functional interpretation of identified genetic variants by NGS platforms is now essential. Available approaches nowadays include bioinformatics, cell and molecular biology and animal models. Recent advances, such as the CRISPR-Cas9, ZFN and TALEN systems, have been already used as a powerful tool with this objective. However, the use of cell lines is of limited value due to the CRC heterogeneity and its close interaction with microenvironment. Access to tridimensional cultures or organoids and xenograft models that mimic the in vivo tissue architecture could revolutionize functional analysis. This review will focus on the application of state-of-the-art functional studies to better tackle new genes involved in germline predisposition to this neoplasm.

Integrated Analysis of Germline and Tumor DNA Identifies New Candidate Genes Involved in Familial Colorectal Cancer.

Colorectal cancer (CRC) shows aggregation in some families but no alterations in the known hereditary CRC genes. We aimed to identify new candidate genes which are potentially involved in germline predisposition to familial CRC. An integrated analysis of germline and tumor whole-exome sequencing data was performed in 18 unrelated CRC families. Deleterious single nucleotide variants (SNV), short insertions and deletions (indels), copy number variants (CNVs) and loss of heterozygosity (LOH) were assessed as candidates for first germline or second somatic hits. Candidate tumor suppressor genes were selected when alterations were detected in both germline and somatic DNA, fulfilling Knudson's two-hit hypothesis. Somatic mutational profiling and signature analysis were also performed. A series of germline-somatic variant pairs were detected. In all cases, the first hit was presented as a rare SNV/indel, whereas the second hit was either a different SNV (3 genes) or LOH affecting the same gene (141 genes). BRCA2, BLM, ERCC2, RECQL, REV3L and RIF1 were among the most promising candidate genes for germline CRC predisposition. The identification of new candidate genes involved in familial CRC could be achieved by our integrated analysis. Further functional studies and replication in additional cohorts are required to confirm the selected candidates.

Mutational Signatures in Cancer (MuSiCa): a web application to implement mutational signatures analysis in cancer samples.

BACKGROUND: Mutational signatures have been proved as a valuable pattern in somatic genomics, mainly regarding cancer, with a potential application as a biomarker in clinical practice. Up to now, several bioinformatic packages to address this topic have been developed in different languages/platforms. MutationalPatterns has arisen as the most efficient tool for the comparison with the signatures currently reported in the Catalogue of Somatic Mutations in Cancer (COSMIC) database. However, the analysis of mutational signatures is nowadays restricted to a small community of bioinformatic experts. RESULTS: In this work we present Mutational Signatures in Cancer (MuSiCa), a new web tool based on MutationalPatterns and built using the Shiny framework in R language. By means of a simple interface suited to non-specialized researchers, it provides a comprehensive analysis of the somatic mutational status of the supplied cancer samples. It permits characterizing the profile and burden of mutations, as well as quantifying COSMIC-reported mutational signatures. It also allows classifying samples according to the above signature contributions. CONCLUSIONS: MuSiCa is a helpful web application to characterize mutational signatures in cancer samples. It is accessible online at http://bioinfo.ciberehd.org/GPtoCRC/en/tools.html and source code is freely available at https://github.com/marcos-diazg/musica .

Zeb1 in Stromal Myofibroblasts Promotes Kras-Driven Development of Pancreatic Cancer.

The transcription factor Zeb1 has been identified as a crucial player in Kras-dependent oncogenesis. In pancreatic ductal adenocarcinoma (PDAC), Zeb1 is highly expressed in myofibroblasts and correlates with poor prognosis. As Kras mutations are key drivers in PDAC, we aimed here to assess the necessity of Zeb1 for Kras-driven PDAC and to define the role of Zeb1-expressing myofibroblasts in PDAC development. Genetically engineered mice with conditional pancreatic KrasG12D and Trp53 mutations (KPC) were crossed with Zeb1 haploinsufficient mice (Z+/-). Extensive PDAC was prominent in all 20-week-old KPC;Z+/+ mice, whereas only low-grade precursor lesions were detected in age-matched KPC;Z+/- littermates, with PDAC developing eventually in KPC;Z+/- aged animals. Zeb1 expression in myofibroblasts occurred early in tumorigenesis and Zeb1 haploinsufficiency retarded native expansion of stromal myofibroblasts during precursor-to-cancer progression. Zeb1 downregulation in mPSC repressed their activated gene profile, impaired their migratory and proliferative activity, and attenuated their tumor-supporting features. Conditioned media from Z+/+ mouse-activated (myofibroblast-like) pancreatic stellate cells (mPSC) boosted Ras activity in pancreatic cancer cells carrying mutant Kras; this effect was not observed when using conditioned media from Z+/- mPSC, revealing a paracrinal cooperative axis between Zeb1-expressing PSC and oncogenic Kras-bearing tumor cells. We conclude that Zeb1-expressing stromal myofibroblasts enable a heterotypic collaboration with the Kras-fated epithelial compartment, thus supporting pancreatic malignancy.Significance: Zeb1 expression in stromal myofibroblasts supports PDAC development via collaboration with the epithelial compartment bearing oncogenic Kras mutations. Cancer Res; 78(10); 2624-37. ©2018 AACR.